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The equilibrium shape of strained InAs quantum dots grown epitaxially on a 
GaAs(OOl) substrate is derived as a function of volume. InAs surface energies are 
calculated within density-functional theory, and a continuum approach is applied 
for the elastic relaxation energies. 



1 Introduction 

The strain-induced self-assembly of small three-dimensional islands during het- 
eroepitaxial growth represents an efficient way of producing quantum dots. 
A frequently studied example is InAs/GaAs(100)(il with the lattice mismatch 
amounting to about 7%. InAs grows on GaAs in the Stranski-Krastanov mode: 
Due to the smaller surface energy of InAs as compared to GaAs first a wet- 
ting layer forms, but when more InAs is deposited coherent, dislocation-free 
three-dimensional islands appear. The driving force for the island-formation 
originates in the gain of elastic relaxation energy which overcompcnsates the 
energetical costs due to the increased surface area. However, the details of the 
growth mechanism, especially the reason behind the narrow siza-distribution 
found for InAs/GaAs(100) quantum dots, are still controversially 

In this paper we will derive the equilibrium shape of the InAs quantum 
dots as a function of size. With restricting ourselves to the energetics we do 
not imply that kinetic effects are unimportant for shaping the dots. However, 
we expect the equilibrium 3D-island shapes to be observable under appropriate 
experimental conditions. When the concentration of dots is low, shape equi- 
libration by atomic diffusion on the small islands will be faster than material 
exchange between the islands leading to Ostwald ripening. 

To calculate the equilibrium 3D-island shapes we have computed InAs sur- 
face energies from first-principles for a variety of surface orientations. A con- 
tinuum approach is adequate to calculate the elastic relaxation energy. Com- 
bining these two contributions we get the total energy, which is minimized 
with respect to the island shape. The more delicate energetical effects, like 
island- island interaction, the strain-dependence of surface energy etc., which 
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have to be considered when discussing island sizesP can be neglected for the 
purpose of predicting island shapes at fixed volume. 

2 InAs Surface Energies 

Surface energies have been calculated using density-functional theory, with the 
local-density approximation being applied to the exchange-correlation func- 
tional. The surface is represented by a periodically repeated slab, about 10 
atomic layers thick, with the topmost 4 layers being fully relaxed. The In and 
As atoms are described by ab initio pseudopotentials, and the electron den- 
sity is calculated using special k-point sets with a density in reciprocal space 
equivalent to 64 k-points in the whole (100) (lxl) surface Brillouin-zone. The 
wave- functions are expanded into plane waves with a kinetic energy fClOJRy. 
We have employed a generalized version of the computer code fhi93cpQ to 
calculate the energy density according to Chetty and MartinB The surface en- 
ergy is derived from integrals of the energy density over appropriately chosen 
subvolumes of the supercell. For a detailed description of similar calculations 
for GaAs see Ref.||. 



Table 1: The equilibrium surface reconstructions of InAs under As-rich conditions and their 

surface energies. 



orientation 


reconstruction 


surface energy [meV/A 2 ] 


(110) 


(lxl) relaxed cleavage plane 


41 


(100) 


a(2x4) 


48 


(100) 


c(4x4) 


47 


(111) 


(2x2) In vacancy 


42 


(ill) 


(2x2) As trimer 


36 



As the InAs surface reconstructions are expected to be similar to those of 
GaAs we have carried out calculations for the same surface reconstructions as 
in Ref. H Epitaxial growth most often takes place under As-rich conditions, 
thus the data in Table [j] refer to surfaces in equilibrium with bulk As ( A7 struc- 
ture). As opposed to GaAs the As-terminated InAs(llO) (lxl) surface does 
not become stable, instead the relaxed cleavage plane is energetically preferred. 
For InAs(100) the As-terminated c(4x4) reconstruction yields the lowest sur- 
face energy. However, the energy difference with respect to the a(2x4) re- 
construction is so small that our calculation is compatible to the experimental 
observation of a (2x4) reconstructionQ Our (111) and (111) equilibrium recon- 
structions are consistent with recent core-level and valence-band photoemission 
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studiesffl The (lll)(2x2) As-trimer reconstruction, which is stable for GaAs 
under As-rich conditions, does not become stable in case of InAs. The surfaces 
in Table [j] are thermodynamically stable against faceting into each other. In 
fact, all four facets have been observed on larg^and thus presumably relaxed, 
InAs islands grown on a GaAs(OOl) substrateEEl 
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Figure 1: Elastic energy per volume E e i a3 t/V vs. surface energy per volume E sur f/V for 
InAs islands with volume V = 2.88 X 10 5 A 3 . Squares: square based pyramid with {101} 
faces and (OOl)-truncated {101}-pyramids. Diamonds: square based pyramids with {111} 
and {111} faces and (OOl)-truncated pyramids. Triangles up: "huts" with {111} and {111} 
faces. Triangles down: square based {101} pyramids with {lll}-truncated edges. Dots: 
islands with {101}, {111}, and {111} faces. Filled symbols denote numerical results, while 
open circles correspond to a simple analytical approximation for (OOl)-truncated "mesa- 
shaped" islands. It is assumed that the elastic energy does not change when the (almost 
fully relaxed) top of an island is cut off. Full lines connect islands that are created in this 
way, varying the height of the (001) surface plane. The dashed line is the curve of constant 
total energy E e i aat + E aur f that selects the equilibrium shape. 



3 Elastic Energy and Equilibrium Shape of Quantum Dots 

The elastic energy is calculated within continuum theory, using a finite-element 
type of approach. The strain field both in the island and in. the substrate is 
fully accounted for. For simplicity, the same elastic constantstil have been taken 
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both for the InAs island and the GaAs substrate. The results are summarized 
in Fig. [I] for a specific volume V of the islands. Generalization to arbitrary 
volume can be done by means of the scaling relations E e i ast ~ V and E sur f ~ 
V 2 / 3 . The equilibrium island-shape is marked by the point where the lines of 
constant total energy E e i ast + E sur j (dashed line in Fig. [j]) touch the manifold 
of island energy-curves from below. The volume-dependence of this shape can 
be inferred from the same figure by rescaling the E sur f/V axis: For larger 
volume V the slope of the total-energy lines becomes smaller. 

The equilibrium island shape results from the competition between elastic 
and surface energy. According to the different scaling properties of E e i ast 
and E sur f for large volumes the elastic energy dominates, favoring a steep 
pyramidal shape, while small islands become more and more flat. Within the 
configuration space we have examined the optimum shapes can be described as 
mesa-type hills bounded by {101}, {111} and {111} faces and an (OOl)-surfape 
on the top. While they differ from the {101}-pyramids seen in experiment 113 
they are in fact similar to shapes observed for InP / GalnP islandslI3 

^From Fig. |l| it is obvious that the shape evolution is continuous with 
respect to the volume. This causes an additional T^-dependence in E e J^ l t {V) 
and E^ l j(V), which does not follow "standard" scaling anymore. Though 
our total-energy expression does not yield any optimum island size by its con- 
struction, this effect is important for a theory of equilibrium island sized 

References 

1. D. Leonard et al, Appl. Phys. Lett. 63, 3203 (1993). 

2. C. Priester and M. Lannoo, Phys. Rev. Lett. 75, 93 (1995). 

3. V. A. Shchukin et al, Phys. Rev. Lett. 75, 2968 (1995). 

4. R. Stumpf and M. Schemer, Comput. Phys. Commun. 79, 447 (1994). 

5. N. Chetty and R. M. Martin, Phys. Rev. B 45, 6074 (1992). 

6. N. Moll et al, Phys. Rev. B. (in press). 

7. H. Yamaguchi and Y. Horikoshi, Phys. Rev. B 51, 9836 (1995). 

8. L. 6. Olsson et al, Phys. Rev. B 53, 4734 (1996). 

9. C. B. M. Andersson et al, Surf. Sci. 347, 199 (1996). 

10. E. Steimetz et al, to be published. 

11. The experimental second and third order elastic moduli of GaAs were 
used. Data in Science and Technology, Semiconductors, Group LV Ele- 
ments and III-V Compounds, Ed. O. Madelung (Springer, Berlin, 1991). 

12. S. Ruvimov et al, Phys. Rev. B 51, 14766 (1995). 

13. K. Georgsson et al, Appl. Phys. Lett. 67, 2981 (1995). 



4 



